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Abstract 

We argue that turbulence in superfluids is governed by two dimen- 
sionless parameters. One of them is the intrinsic parameter q which 
characterizes the friction forces acting on a vortex moving with respect 
to the heat bath, with q~ l playing the same role as the Reynolds num- 
ber Re = UR/v in classical hydrodynamics. It marks the transition 
between the " laminar" and turbulent regimes of vortex dynamics. The 
developed turbulence described by Kolmogorov cascade occurs when 
Re S> 1 in classical hydrodynamics, and q <C 1 in the superfluid hy- 
drodynamics. Another parameter of the superfluid turbulence is the 
superfluid Reynolds number Re s = UR/k, which contains the cir- 
culation quantum k characterizing quantized vorticity in superfluids. 
This parameter may regulate the crossover or transition between two 
classes of superfluid turbulence: (i) the classical regime of Kolmogorov 
cascade where vortices are locally polarized and the quantization of 
vorticity is not important; and (ii) the quantum Vinen turbulence 
whose properties are determined by the quantization of vorticity. The 
phase diagram of the dynamical vortex states is suggested. 

PACS numbers: 67.40.Vs, 43.37.+q, 4732. Cc, 67.57.Fg 
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1 Introduction 



The hydrodynamics of superfluid liquid exhibits new features with respect 
to conventional classical hydrodynamics, which are important when the tur- 
bulence in superfluids is considered [1]. 

(i) The superfluid liquid consists of two mutually penetrating components 
- the frictionless superfluid component and the viscous normal component. 
That is why different types of turbulent motion are possible depending on 
whether the normal and superfluid components move together or separately. 
Here we are interested in the most simple case when the dynamics of the 
normal component can be neglected. This occurs, for example, in superfluid 
phases of 3 He where the normal component is so viscous that it is practically 
clamped by the container walls. Its role is to provide the preferred reference 
frame, where the normal component and thus the heat bath are at rest. The 
turbulence in the superfluid component with the normal component at rest 
is referred to as the superfluid turbulence. 

(ii) The important feature of the superfluid turbulence is that the vor- 
ticity of the superfluid component is quantized in terms of the elementary 
circulation quantum k. So the superfluid turbulence is the chaotic motion 
of well determined and well separated vortex filaments [1]. Using this we 
can simulate the main ingredient of the classical turbulence - the chaotic 
dynamics of the vortex degrees of freedom of the liquid. 

(iii) The further simplification comes from the fact that the dissipation 
of the vortex motion is not due to the viscosity term in the Navier-Stokes 
equation which is proportional to the velocity gradients V 2 v in classical 
liquid, but due to the friction force acting on the vortex when it moves with 
respect to the heat bath (the normal component). This force is proportional 
to velocity of the vortex, and thus the complications resulting from the V 2 v 
term are avoided. 

Here we discuss how these new features could influence the superfluid 
turbulence. 

2 Coarse-grained hydrodynamic equation 

The coarse-grained hydrodynamic equation for the superfluid vorticity is 
obtained from the Euler equation for the superfluid velocity v = v s after 
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averaging over the vortex lines (see the review paper [2]). Instead of the 
Navier-Stokes equation with V 2 v term one has 

_ + V/i = v x uj - (1) 

— a'(v — v n ) x uj + a uj x (uj x (v - v n )) . (2) 

Here v n is the velocity of the normal component; uj = V x v is the superfluid 
vorticity; uj = uj/uj; and dimensionless parameters a' and a come from the 
reactive and dissipative forces acting on a vortex when it moves with respect 
to the normal component. These parameters are very similar to the Hall 
resistivity p xy and p xx in the Hall effect. For vortices in fermionic systems 
(superfluid 3 He and superconductors) they were calculated by Kopnin [3], 
and measured in 3 He-B in the broad temperature range by Bevan et. al. 
[4] (see also [5], where these parameters are discussed in terms of the chiral 
anomaly) . 

The terms in Eq.(l) are invariant with respect to the transformation 
v — > v(r — ut) + u as in classical hydrodynamics. However, the terms in 
Eq.(2) are not invariant under this transformation, since there is the preferred 
reference frame in which the normal component is at rest. They are invariant 
under the full Galilean transformation when the normal component is also 
involved: v — > v(r — ut) + u and v n — > v n + u. 

Further we shall work in the frame where v n = 0, but we must remember 
that this frame is unique. In this frame the equation for superfluid hydrody- 
namics is simplified: 

<9v _ 

— + V/i = (1 — a')v xw + aux(wxv). (3) 

C6 

After rescaling the time, t — (1 — a')t, one obtains equation which depends 
on a single parameter q — a/(l — a'): 

+ V/i = v x uj + q uj x (uj x v) . (4) 

Now the first three terms together are the same as inertial terms in classical 
hydrodynamics. They satisfy the modified Galilean invariance (in fact the 
transformation below changes the chemical potential, but this does not influ- 
ence the vortex degrees of freedom which are important for the phenomenon 
of turbulence) : 

v(t, r) -> v(t, r - ut) + u . (5) 
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On the contrary the dissipative last term with the factor q in Eq.(4) is not 
invariant under this transformation. This is in contrast to the conventional 
liquid where the whole Navier-Stokes equation which contains viscosity 

— +V/i = vxw + z/V 2 v , (6) 
at 

is Galilean invariant, and where there is no preferred reference frame. 

Such a difference between the dissipative last terms in Eqs. (6) and (4) 
is very important: 

(1) The role of the Reynolds number, which characterizes the ratio of 
inertial and dissipative terms in hydrodynamic equations, in the superfluid 
turbulence is played by the intrinsic parameter 1/q. This parameter does 
not depend on the characteristic velocity U and size R of the large-scale flow 
as distinct from the conventional Reynolds number Re = RU jv in classical 
viscous hydrodynamics. That is why the turbulent regime occurs only at 
1/q > 1 even if vortices are injected to the superfluid which moves with large 
velocity U. This rather unexpected result was obtained in recent experiments 
with superfluid 3 He-B [6]. 

(2) In the conventional turbulence the large-scale velocity U is always 
understood as the largest characteristic velocity difference in the inhomo- 
geneous flow of classical liquid [7]. In the two-fluid system the velocity U 
is the large-scale velocity of superfluid component with respect to the nor- 
mal component, and this velocity (the so-called counterflow velocity) can be 
completely homogeneous. 

(3) As a result, as distinct from the classical hydrodynamics, the energy 
dissipation which is produced by the last term in Eq.(4) depends explicitly 
on U : 

e = —E = -(v • —=) = -q{v ■ (<2> x (u x v))) ~ quU 2 . (7) 
at 

(4) The onset of the superfluid turbulence was studied in Ref. [8], where 
the model was developed which demonstrated that the initial avalanche- 
like multiplication of vortices leading to turbulence occurs when q < 1 in 
agreement with experiment [6]. The existence of two regimes in the initial 
development of vorticity is also supported by earlier simulations by Schwarz 
who noted that when a (or q) is decreased the crossover from a regime of 
isolated phase slips to a phase-slip cascades and then to the fully developed 
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vortex turbulence occurs [9]. One can expect that the well developed turbu- 
lence occurs when g < 1, and here we shall discuss this extreme limit. In 
3 He-B the condition q <C 1 can be ralized at temperatures well below 0.6T C 
[6]. However, we do not consider a very low T where instead of the mutual 
friction the other mechanisms of dissipation take place such as excitation of 
Kelvin waves [10] and vortex reconnection [11]. The latter leads to formation 
of cusps and kinks on the vortex filaments whose fast dynamics creates the 
burst of different types of excitations in quantum liquids: phonons, rotons, 
Kelvin waves and fermionic quasiparticles. The burst of gravitational waves 
from cusps and kinks of cosmic strings was theoretically investigated by the 
cosmological community (see e.g. [12]), and the obtained results are very 
important for the superfluid turbulence at a very low temperature. 

(5) We expect that even at q <C 1 two different states of turbulence are 
possible, with the crossover (or transition) between them being determined 
by q and by another dimensionless parameter Re s = UR/k, where k is the 
circulation around the quantum vortex. The coarse-grained hydrodynamic 
equation (4) is in fact valid only in the limit Re s ^> 1, since the latter 
means that the characteristic circulation of the velocity r = UR of the large- 
scale flow substantially exceeds the circulation quantum k, and thus there 
are many vortices in the turbulent flow. When Re s decreases the quantum 
nature of vortices becomes more pronounced, and we proceed from the type 
of the classical turbulence which is probably described by the Komogorov 
cascade, to the quantum regime which is probably described by the Vinen 
equations for the average vortex dynamics [13]. 

Let us consider the possibility of the Kolmogorov state of the superfluid 
turbulence. 

3 Kolmogorov cascade 

In classical turbulence, the large Reynolds number Re — UR/v >> 1 leads to 
the well separated length scales or wave numbers. As a result the Kolmogorov- 
Richardson cascade takes place in which the energy flows from small wave 
numbers k m i n ~ 1/R (large rings of size R of the container) to high wave 
number k = l/r where the dissipation occurs. In the same manner in our 
case of the superfluid turbulence the necessary condition for the Kolmogorov 
cascade is the big ratio of the inertial and dissipative terms in Eq.(4), i.e. 
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l/g» 1. 

In the Kolmogorov-Richardson cascade, at arbitrary length scale r the 
energy transfer rate to the smaller scale, say r/2, is e = E r /t r , where E r = v 2 
is the kinetic energy at this scale, and t r = r/v r is the characteristic time. 
The energy transfer from scale to scale must be the same for all scales, as a 
result one has 

E r vl U 3 
e = — = — = constant = — — . (8) 
t r r R y J 

From this equation it follows that 

v r = e 1 / 3 r 1 / 3 . (9) 

This must be valid both in classical and superfluid liquids [14]. What is 
different is the parameter e: it is determined by the dissipation mechanism 
which is different in two liquids. 

From Eq.(7) with u r = v r /r it follows that as in the classical turbulence 
the main dissipation occurs at the smallest possible scales, but the structure 
of e is now different. Instead of e = z/tVo/ r o 1X1 classical liquids, we have now 

e ~ qu r0 U 2 ~ qU 2 ^ = qU 2 e 1 ' 3 r 2/3 . (10) 

Since e = U 3 /R one obtains from Eq.(10) that the scale ro at which the main 
dissipation occurs and the characteristic velocity v ro at this scale are 

r^q 3 ' 2 R , v r0 ~q 1/2 U. (11) 

This consideration is valid when r R and v r0 U, which means that 
1/q ^> 1 is the condition for the Kolmogorov cascade. In classical liquids 
the corresponding condition for the well developed turbulence is Re ^> 1. In 
both cases these conditions ensure that the kinetic terms in the hydrodynamic 
equations are much larger than the dissipative terms. In the same manner 
as in classical liquids the condition for the stability of the turbulent flow is 
Re > 1, one may suggest that the condition for the stability of the discussed 
turbulent flow is 1/q > 1. This is supported by observations in 3 He-B where 
it was demonstrated that at high velocity U but at q > 1 the turbulence is 
not developed even after vortices were introduced into the flow [6]. 
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As in the Kolmogorov cascade for the classical liquid, in the Kolmogorov 
cascade of superfluid turbulence the dissipation is concentrated at small 
scales, 

e~qU 2 [ R ^ V -L~ qU ^, (12) 
Jr r r r 

while the kinetic energy is concentrated at large scale of container size: 

E = [ R -v 2 = f R -(erf/ 3 = (ei?) 2 / 3 = U 2 . (13) 

Jr T Jr T 

The dispersion of the turbulent energy in the momentum space is the 
same as in classical liquid 



l- (er)/ =l TW^L ikE{k) < 

E(k) = e" 3 k-*' 3 . (14) 



As distinct from the classical liquid where ko is determined by viscosity, 
in the superfluid turbulence the cut-off ko is determined by mutual friction 
papameter q: k — l/r — R~ l q~^l 2 . 



4 Crossover to Vinen quantum turbulence 

At a very small q the quantization of circulation becomes important. The 
condition of the above consideration is that the relevant circulation can be 
considered as continuous, i.e. the circulation at the scale tq is larger than 
the circulation quantum: f r o?"o > K - This gives 

UR 

Vro r o — q 2 UR = q 2 nRe s > k , Re s = , (15) 

K 

i.e. the constraint for the application of the Kolmogorov cascade is 

Re s > \ > 1 . (16) 
q 2 

Another requirement is that the characteristic scale r must be much larger 
that the intervortex distance /. The latter is obtained from the vortex density 
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in the Kolmogorov state % = l~ 2 = uJto/k = ^r-o/( r o K )- The condition I <C tq 
leads again to the equation t> r o r o > K an d thus to the criterion (16). 

Note that here for the first time the 'superfluid Reynolds number' Re s 
appeared, which contains the circulation quantum. Thus the superfluid 
Reynolds number is responsible for the crossover or transition from the clas- 
sical superfluid turbulence, where the quantized vortices are locally aligned 
(polarized), and thus the quantization is not important, to the quantum 
turbulence developed by Vinen. 

We can now consider the approach to the crossover from the quantum 
regime - the Vinen state which probably occurs when Re s g 2 < 1. Accord- 
ing to Vinen [13] the characteristic length scale, the distance between the 
vortices or the size of the characteristic vortex loops, is determined by the 
circulation quantum and the counterflow velocity, / = Xk/U, where A is the 
dimensionless intrinsic parameter, which probably contains a' and a. The 
vortex density in the Vinen state is 

nv = I' 2 - A 2 ^ = ^Re 2 , (17) 



K 



It differs from the vortex density in the Kolmogorov state 

VrQ U 1 Re s 
nK = ^ ~ q^R = ' (18) 

which depends not only on the counterflow velocity U, but also on the con- 
tainer size R. 

If the crossover between the classical and quantum regimes of the tur- 
bulent states occurs at Re s g 2 = 1, the two equations (17) and (18) must 
match each other in the crossover region. But this occurs only if A 2 = q. If 
A 2 7^ q there is a mismatch, and one may expect that either the two states 
are separated by the first-order phase transition, or there is an intermediate 
region where the superfluid turbulence is described by two different micro- 
scopic scales such as r and I. Based on the above consideration one may 
suggest the following phase diagram of different regimes of collective vortex 
dynamics in Fig. 1. 
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Figure 1: Possible phase diagram of dynamical vortex states in Re s , q plane. 
At large flow velocity Re s 3> 1 the boundary between the turbulent and 
'laminar' vortex flow approaches the vertical axis q = go ~ 1 as suggested 
by experiment [6]. The thick line separates the developed turbulence of the 
classical type, which is characterized by the Kolmogorov cascade at small q, 
and the quantum turbulence of the Vinen type. 
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5 Discussion 



The superfluid turbulence is the collective many-vortex phenomenon which 
can exist in different states. Each of the vortex states can be characterized by 
its own correlation functions. For example, the states can be characterized 
by the loop function 



In the limit when the length L of the loop C is much larger than the inter- 
vortex distance I one may expect the general behavior g(L) ~ exp(— (L/7) 7 ) 
where the exponent 7 is different for different vortex states. 

One can expect the phase transitions between different states of of col- 
lective vortex dynamics. One of such transitions which appeared to be 
rather sharp has been observed between the 'laminar' and 'turbulent' dy- 
namics of vortices in superfluid 3 He-B [6]. It was found that such transi- 
tion was regulated by intrinsic velocity-independent dimensionless parameter 
q — a/(l — a'). However, it is not excluded that both dimensionless param- 
eters a and a' are important, and also it is possible that only the initial 
stage of the formation of the turbulent state is governed by these parame- 
ters [8]. Another transition (or maybe crossover) is suggested here between 
the quantum and classical regimes of the developed superfluid turbulence, 
though there are arguments that the classical regime can never be reached 
because the vortex stretching is missing in the superfluid turbulence [15]. 

In principle the parameters a and a' may depend on the type of the 
dynamical state, since they are obtained by averaging of the forces acting on 
individual vortices. The renormalization of these parameters a(L) and cx'(L) 
when the length scale L is increasing may also play an important role in the 
identification of the turbulent states, as in the case of the renormalization- 
group flow of similar parameters in the quantum Hall effect (see e.g. [16]). 

I thank V.B. Eltsov, D. Kivotides, N.B. Kopnin and M. Krusius for dis- 
cussions. This work was supported by ESF COSLAB Programme and by the 
Russian Foundations for Fundamental Research. 
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